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Multiferroics permit the magnetic control of 
the electric polarization and electric control of 
the magnetization 1 8 . These static magnetoelec¬ 
tric (ME) effects are of enormous interest: The 
ability to read and write a magnetic state current- 
free by an electric voltage would provide a huge 
technological advantage 5 8 . Dynamic or optical 
ME effects are equally interesting because they 
give rise to unidirectional light propagation as 
recently observed in low-temperature multifer- 
roics^H^. This phenomenon, if realized at room 
temperature, would allow the development of op¬ 
tical diodes which transmit unpolarized light in 
one, but not in the opposite direction. Here, we 
report strong unidirectional transmission in the 
room-temperature multiferroic BiFeC>3 over the 
gigahertz—terahertz frequency range. Support¬ 
ing theory attributes the observed unidirectional 
transmission to the spin-current driven dynamic 
ME effect. These findings are an important step 
toward the realization of optical diodes, supple¬ 
mented by the ability to switch the transmission 
direction with a magnetic or electric field. 

BiFeOs is by the far most studied compound among 
multiferroic and magnetoelectric materials. While ex¬ 
perimental studies have already reported about the first 
realizations of the ME memory function using BiFeOs 
based device^®, the origin of the ME effect is still under 
debate due to the complexity of the material. Because 
of the low symmetry of iron sites and iron-iron bonds, 
the magnetic ordering can induce local polarization via 
each of the three canonical term ^ 5 - the spin-current, 
exchange-striction and single-ion mechanisms. While the 
spin-current term has been identified as the leading con¬ 
tribution to the magnetic ally i nduced ferroelectric polar¬ 
ization in various studies 4 the spin-driven atomic 

displacements 18 and the electrically induced shift of the 
spin-wave (magnon) resonance! 19 ! were interpreted based 
on the exchange-striction and single-ion mechanisms, re¬ 
spectively. 

In the magnetically ordered phase below T/v=640 K, 
BiFeOs possesses an exceptionally large spin-driven po¬ 
larization 1 ^, if not the largest among all known multi¬ 
ferroic materials. Nevertheless, its systematic study has 
long been hindered by the huge lattice ferroelectric po¬ 
larization (Pq) developing along one of the cubic (111) 
directions at the Curie temperature 7^=1100K and by 


the lack of single-domain ferroelectric crystals. Owing 
to the coupling between P 0 and the spin-driven polar¬ 
ization, in zero magnetic field they both point along the 
same [ 111 ] axis. A recent systematic study of the static 
ME effect revealed additional spin-driven polarization or¬ 
thogonal to [ilijni. 

The optical ME effect of the magnon modes in multifer¬ 
roics, which gives rise to the unidirectional transmission 
in the gigahertz-terahertz frequency range, has recently 
become a hot topic in materials science. The difference in 
the absorption coefficients (a) of beams counterpropagat- 
ing in such ME media—called directional dic hroi sm—can 
be expressed for linear light polarization as ® 0 

A a k (u) = a +fe (w) - a_ fc (w) « - Xs™ Ml- 

( 1 ) 

The dynamic ME susceptibility tensors x me (u;) and 
X ern (uj ) respectively describe the magnetization gener¬ 
ated by the oscillating electric field of light, A= 
(5o//io) 1 ^ 2 X^ e ( CJ )F^, and the electric polarization in¬ 
duced by its oscillating magnetic field, A P$ = 
(eoMo) 1 / 2 X^ Here Sq and /xo are the vacuum 

permittivity and permeability, respectively, while 7 and S 
stand for the Cartesian coordinates. Since the two cross- 
coupling tensors are connected by the time-reversal op¬ 
eration [...]' according to [x^s^)]' = — the di¬ 
rectional dichroism becomes Aak(co) = — 

[x^f( ( ^)] / }- I n other words, the directional dichroism 
emerges for simultaneously electric- and magnetic-dipole 
active excitations and its magnitude is determined by the 
time-reversal odd parts of the off-diagonal X^ e (^) t ensor 
elements 10 11 2 Hf 22 . The schematic representation of the 
optical diode function in ME media is shown in Fig. []] 

In the cycloidal spin state of BiFeOs ^1, several low- 
frequency collective modes have been observed by spec¬ 
troscopic methods including light absorptiorP ^ 25 and 
Raman spectroscopy 6 19 26 . Though the electric-field- 
induced shift of the resonance frequencies observed in the 
Raman study indicates the ME nature of these magnon 
mode!®, the optical ME effect has not been investi¬ 
gated in BiFeOs. Here, we performed absorption mea¬ 
surements in the gigahertz-terahertz spectral range on 
single-domain ferroelectric BiFeOs crystal^! with Po 
along [111] between room temperature and T=4K in 
magnetic fields up to juoff=17T. We found that some of 
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FIG. 1: | Optical diode function in multiferroics. 

Ferro-type ordering of the local electric dipoles (red arrows) 
and magnetic moments (green arrows) produces a ferroelec¬ 
tric polarization P and a spontaneous magnetization M, re¬ 
spectively. Light interacts with both ferroic order parameters, 
hence, upon illumination P and M oscillate coherently with 
the electromagnetic field around their equilibria. The light in¬ 
duced polarization has contributions from both the usual di¬ 
electric permittivity and the optical ME effect x em (^)- While 
the first contribution is independent of the light propagation 
direction, the polarization induced via the optical ME effect 
has opposite sign for counter-propagating light beams. This 
can give rise to either a constructive or a destructive inter¬ 
ference between the two terms. Similarly, the magnetization 
dynamics is governed by the interference between the magne¬ 
tization induced via the magnetic permeability and the optical 
ME effect x me (^)- Consequently, the transmitted intensity 
depends on the propagation direction (intense and pale yel¬ 
low beams) even for unpolarized light and can be exploited to 
produce optical diodes transmitting light in one, but not in 
the opposite direction. The transmitting direction can be re¬ 
versed by switching the sign of either P via an electric voltage 
(V) or M by an external magnetic field (H). 


the magnon modes exhibit strong unidirectional trans¬ 
mission. We identified the minimal set of spin-driven- 
polarization terms and quantitatively reproduced both 
the spectral shape and the field dependence of the direc¬ 
tional dichroism solely by the spin-current mechanism. 

The experimental configurations are schematically il¬ 
lustrated in Fig. [2] Absorption spectra were obtained for 
light beams propagating along [001] with two orthogonal 
linear polarizations, E w ||[110] and E w ||[110]. Static mag¬ 
netic fields (=b H) were applied perpendicular to the light 
propagation direction along either [110] or [110]. 

In simple magnets, such as ferromagnets, the sign 
change of the magnetization corresponds to the time re¬ 
versal operation. Thus, it is equivalent to the reversal of 
the light propagation direction. Owing to experimental 
limitations, in such cases, the absorption change upon 



FIG. 2: | Experimental configurations used to detect 

unidirectional transmission in BiFeOs. a, Pseudocu- 
bic unit cell of BiFeOs showing the positions of Bi, Fe and 
O ions. The lattice ferroelectric polarization, Po||[111], is 
schematically indicated on the Fe site, b, Illustration of the 
three equivalent directions of the cycloidal ordering vector q i 
on the Fe sublattice. The frame of reference is common to 
all panels, c, In magnetic fields (±1L) applied along [110], 
cycloidal domains with c \2 and q 3 are e qually favoured, while 
the domain with qi is suppresse d 28 * 29 *, d, In magnetic fields 
(±1L) applied along [110], only the cycloidal domain with qi 
is stabkP^^. e & f , The propagation direction (k) and the 
two orthogonal polarizations of light beams traveling in the 
material. 


the magnetic field induced reversal of the magnetization, 
A OLH=OL+H,+k-ot—H,+k, is typically detected instead of 
the absorption change associated with the reversal of 
the light propagation direction, A ak=C£+H,+k—ct+H,-k- 
Though the relation Aa^ = A<a# does not necessar¬ 
ily hold for complex spin structures, such as BiFeOs, 
Act'/ C and Aa# spectra obtained from our calculations 
are equal within 1-2 % for the experimental configura¬ 
tions studied here. 

Figure [3] shows the absorption spectra measured in 
four different configurations, i.e. for two orientations 
of the magnetic field and two light polarizations. The 
absorption coefficient at several magnon resonances de¬ 
pends on the sign of the magnetic field. This differ¬ 
ence is stronger for H|| [110] and most pronounced for the 
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FIG. 3: | Absorption spectra of BiFeC >3 in the range of magnon resonances, a-d Magnetic field dependent part 

of the absorption spectra measured at T = 2.5 K in four different configurations, i.e. for the two orientations of the magnetic 
field (H) and the two orthogonal polarizations schematically shown in Fig. [5] The light propagation direction is common 
to all experimental configurations, k|| [001] . Absorption spectra measured in different magnetic fields are shifted vertically in 
proportion to the magnitude of the field, and spectra recorded in +H and —H are plotted with red and black lines, respectively. 
Spectra shown in panel a & c represent absorption form the qi cycloidal domain stabilized by H||[110], while spectra in panel 
b & d have contributions from q 2 and q 3 domains favoured by H || [110]. 




lowest-frequency mode \ko in Fig. [ 3^1 when H|| [110] and 
|| [110]. With increasing magnetic field this resonance 
becomes almost transparent for +iif, while its absorption 
increases for —H. We also measured the absorption spec¬ 
tra with both light polarizations for H|| [001] and could 
not detect any difference between AH. 

In order to reproduce the observed directional dichro- 
ism on a microscopic basis, we adopt the spin model 
of Refs. |28l29] . which successfully describes the mag¬ 
netic field dependence of the magnon resonances (see the 
Methods section). Similarly to the static ME effect, all 
the three basic mechanisms—the spin-current, exchange- 
striction and single-ion mechanism—can in principle con¬ 
tribute to the optical ME effect. By including all 
symmetry-allowed spin-driven polarization terms, we 
calculated the optical ME susceptibilities x me (u) and 


X em (u ;), the diele ctric p ermittivity i(cu) and the magnetic 
permeability fi(cj) 14 20 . Next, we numerically solved the 
Maxwell equations by including these response functions 
in the constitutive relations and calculated the transmis¬ 
sion of linearly polarized incoming beams for both back¬ 
ward and forward propagation. The same calculation 
was done for both field directions, db H. As already men¬ 
tioned, we found Ao^ « Aan irrespective of the magni¬ 
tude of H. 


To identify the spin-driven polarization terms relevant 
to the optical ME effect, we performed a systematic fit¬ 
ting of the measured A qlh{w) by treating the magnitude 
of the different terms as free parameters. We found that 
the directional dichroism spectra are closely reproduced 






















































































































4 


by the following two types of spin-current terms 

P » C = y E{4 1) [e iJ x(S i xS i )] a +(-l)^Ai 2 )[S i xS i ] Q } 
(ij) 

( 2 ) 

where the summation goes over neighbouring spins con¬ 
nected by unit vectors and the integer rii labels 
the hexagonal layers along [111]. The dynamic ME ef¬ 
fect generated by the spin-current terms is described by 
the coupling constants A^ and Xa \ where a=x',y',z' 
stands for the three coordinates along the axes x'Hq*, 
y'lKPoXq*) and z'||P 0 (see Fig.^). 

Figure [4] shows the comparison between the mea¬ 
sured and calculated directional dichroism spectra 
for H || [110] with the two orthogonal light po¬ 
larizations, || [110] and E w ||[110]. The best 
fit was obtained with three independent param¬ 
eters: A^V=0, A^=—2A^1^57.0±3.1 nC/cm 2 , 

=A «34.5±2.4 nC / cm 2 , A^ «11.8±2.9 nC / cm 2 . 
The population of the two cycloidal domains with q 2 
and q3 propagation vectors was kept equal 28 29 . We 
note that this limited set of parameters provides only 
a semi-quantitative description of the mean absorption 
spectra, a(u)=a+ H ^ k (uj)+a- H ^ k (uj). 

We found that additional terms did not further im¬ 
prove the quality of the fit. Hence, the optical ME 
effect in BiFeOs is dominated by two types of spin- 
current polarizations, while the exchange-striction and 
single-ion polarization terms do not significantly con¬ 
tribute to it. This stems from the general nature of 
the spin dynamics in BiFeOs. Due to the very weak 
on-site anisotropy acting on the S = 5/2 iron spins, 
each magnon mode corresponds to pure precessions of 
the spins, where the oscillating component of the spin 
on site i, is perpendicular to its equilibrium direc¬ 

tion, S 9 . This is in contrast to the spin stretching modes 
observed in highly anisotropic magnets 21 . Since neigh¬ 
bouring spins are nearly collinear in the cycloidal state 
with extremely long (62 nm) pitctP^, a dynamic polar¬ 
ization is efficiently induced via spin-current terms such 
as <5P“ oc S? x £S[/ +1 . In contrast, the dynamic po¬ 
larization generated by exchange-striction terms such as 
SP^ oc S 9 • £S[/ +1 is nearly zero. 

Despite its success in quantitatively describing the di¬ 
rectional dichroism spectra observed for H || [110], our 
model may not be complete. When light propagates 
along [001], we predict that directional dichroism should 
be absent for a magnetic field along [ 7777 ft]. While this 
is in agreement with A<a#=0 found for H|| [001], it can¬ 
not account for the finite directional dichroism discerned 
in Figs. [3^ and |3]i for H 11 [110]. This discrepancy may 
come from additional anisotropy terms, neglected in the 
microscopic spin Hamiltonian adopted from Refs. [28129] . 
which further reduce the symmetry of the magnetic state 
and allow the weak directional dichroism observed for 
H|| [110]. 

Finally, we turn to the temperature dependence of the 
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FIG. 4: | Directional dichroism spectra of BiFeC >3 in 

the range of magnon resonances, a &; b, Magnetic field 
dependence of the directional dichroism spectra measured as 
Aq!h(w) at T— 2.5K with the two orthogonal polarizations 
E^||[110] and E w ||[110], respectively. Spectra obtained in dif¬ 
ferent magnetic fields are shifted vertically in proportion to 
the magnitude of the field, which was applied along [110]. 
The field values (common to each panel) are indicated with 
labels on the top of the spectra in panel a. c & d, Directional 
dichroism spectra predicted by our model for the case of pan¬ 
els a & b, respectively. The calculated mode frequencies are 
indicated by dashed lines. For the assignment of the different 
modes see Refs. [28129] . 

directional dichroism presented in Fig.JHJfor HHE^H[110]. 
With increasing temperature the magnon modes soften 25 
and both the mean absorption and the directional dichro¬ 
ism are reduced. Nevertheless, the modes and 

<£>2 ’ ; still exhibit considerable directional dichroism, 
Zk<a#~ 5 cm -1 a t room temperature. At low tempera¬ 
tures, almost perfect unidirectional transmission was ob¬ 
served for the lowest-energy mode T 0 with orthogonal 
light polarization (E w || [110]). Though we expect the 
same at room temperature, Tq is out of our limited spec¬ 
tral window at high temperatures. 

The emergence of strong optical ME effect and the 
corresponding unidirectional transmission require the si¬ 
multaneous breaking of the space- and time-inversion 
symmetries by the coexistence of ferroelectricity and 
(anti)ferromagnetism. While these optical ME phenom¬ 
ena have been investigated recently in various materi¬ 
als hosting multiferroicity at low temperatures, here we 
studied the unidirectional transmission in the spin ex¬ 
citation spectrum of BiFeOs^®, the unique multiferroic 
compound offering a real potential for room tempera¬ 
ture applications up to date. We found that the optical 
ME effect in BiFeOs is robust enough to generate con- 
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Frequency (THz) 



FIG. 5: | Temperature dependence of the directional 

dichroism in BiFeC> 3 . Directional dichroism spectra mea¬ 
sured in /xo!7=±12T at T— 4, 150 and 300 K. The magnetic 
field was applied along [110] and E w ||[110]. Modes and 
<I >2 ’ ; soften and get weaker with increasing temperature, but 
are still clearly observable even at 300 K. The \Eq is visible 
until 150 K while the cannot be detected reliably already 
at 150 K. 

siderable directional dichroism in the gigahertz-terahertz 
range even at room temperature. Based on the current 
progress achieved in the electric control of the magneti¬ 
zation in BiFeOs, we expect that the magnetic switching 
of the transmission direction, demonstrated here, can be 
complemented by the electric control of the optical ME 
effect. Because these functionalities exist at room tem¬ 
perature, they can pave the way for the development of 
optical diodes with electric and/or magnetic control. 


Methods 

Absorption measurements in the terahertz fre¬ 
quency range. The terahertz spectroscopy system con¬ 
sists of a Martin-Puplett type interferometer with a Si 
bolometer operating at 300 mK and a mercury lamp. At 
high temperatures the spectral window of the measure¬ 
ment was limited by the strong radiation load on the 
detector. The light is directed to the sample using light 
pipes. The sample is located in the He exchange gas filled 
sample chamber, which is placed into the cold bore of a 
17 T superconducting solenoid. 

The measurement sequence was started by applying 
high magnetic fields (>12T) at 4K for tens of minutes. 
For H || [110] and [110], this procedure respectively pop¬ 
ulates a single magnetic domain with qi and two domains 
with q 2 and q 3 . Next, spectra were measured in different 
±i7 fields. We did not see any change in the magnetic do¬ 
main population when the —17 T field was applied after 
+17T. 

The zero field absorption spectrum was subtracted 
from the spectra measured in finite fields. This pro¬ 


cedure cancels out diffraction and interference effects 
caused by the sample. The differential absorption coeffi¬ 
cient a(H) — a(0) = — ln[7(i7)/7(0)]d _1 , where 7(0) and 
7(i7) are light intensity spectra in zero and H field. The 
lower envelope of the whole set of differential absorption 
spectra measured in different fields was used to calcu¬ 
late the zero field spectrum. Magnetic field dependent 
absorbtion spectra were evaluated as a sum of the zero 
field spectrum and the corresponding differential spec¬ 
tra. Note that by this method only the field-dependent 
part of this absorption is recovered and field-independent 
features are not captured. While this can cause an ambi¬ 
guity of the mean absorption, the directional dichroism 
spectrum Aan = a(-\-H) — a(—H) are free of such un¬ 
certainties. 

Theoretical calculations. The cycloid of BiFeOs is 
controlled by two Dzyaloshinskii-Moriya (DM) interac¬ 
tions and an easy-axis anisotropy K along the ferroelec¬ 
tric polarization Pq. Whereas the DM interaction D\ 
perpendicular to Po is responsible for the formation of 
the long-period (62 nm) cycloid^ED the DM interaction 
D2 along Po is responsible for a small cycloidal tilt 31 . 

In a magnetic field H, the spin state and the magnon 
excitations of BiFeOs were evaluated from the micro¬ 
scopic Hamiltonian 

n = -JiJ2 s i- s i - J 2 F s* • Si 

(hi) (ijy 

+D\ ^(z' x e *.j) • ( S * x S i) 

(hi) 

+ d 2 • (®* X Si) 

(hi) 

-^^(z'-S i ) 2 -2 MB H-^S i , (3) 

i i 

Ji and J 2 are the nearest and the second nearest neigh¬ 
bour interactions, respectively. The D 1 sum, first pro¬ 
posed by Katsura and coworkers 1 ^, is uniform over the 
lattice, while the sum alternates sign from one hexag¬ 
onal layer to the next. 

The nearest- and next-nearest neighbor exchange inter¬ 
actions Ji=-5.32meV and J2=-0.24meV were obtained 
from recent inelastic neutron scattering measurements 32 . 
Di«0.18meV is determined from the cycloidal pitch. 
D 2 =0.085 meV and 77=0.0051 meV were obtained by fit¬ 
ting the four magnon modes observed in zero fiel d 24 * 25 *. 
The spin state of BiFeOs is solved by using a trial spin 
state that contains harmonics of the fundamental order¬ 
ing wavevector Q. We then minimize the energy (H) 
over the variational parameters of that state. A 1/S 
expansion about the classical limit is used to evaluate 
the mag non mode frequencies at Q as a function of 
fiel d 28 * 29 * . We calculated the optical ME susceptibilities, 
X me (cj) and x em (cj), the dielectric permittivity, £( cj), and 
the magn etic permeability, /7 (c<j) , using the Kubo for- 
mulsP^o] 

To determine the directional dichroism, the Maxwell 
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equations were numerically solved for linearly polarized 
monochromatic plane waves with ±k. Polarization rota¬ 
tion of the transmitted beam was found negligible. We 
also evaluated the directional dichroism spectra using the 
approximate formula in Eq. (1), which is valid when the 
polarization rotation of light can be neglected^. Eq. (1) 
and the numerical solution of the Maxwell equations pro¬ 
vide nearly equivalent Aak spectra if parameters x me (u;), 
X em (cj), i(u) and ju(uj) realistic to BiFeOs are used. Fits 
to the directional dichroism spectra indicate that the po¬ 
larization induced by the spin current associated with D i 
and D 2 can be written in the form given by Eq.(2), where 
and are the dynamic ME couplings for the two 
types of terms, respectively. 
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